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Abstract. Depletion of the liquid density near a solid surface with a weak long- 
ed . range fluid-surface interaction was studied by computer simulations of the liquid- 
\ vapor coexistence of a LJ fluid confined in slitlike pores. In a wide temperature 
range the liquid density decreases towards the surface without the formation of a 
vapor layer between the liquid and the solid surface. This evidences the absence of a 
I ' drying transition up to the liquid-vapor critical point. Two contributions to the excess 

desorption Y to t were found. The first one ~ pbuik £ exists at any temperature and 
diverges as the bulk correlation length £ when approaching the liquid-vapor critical 
^ ' temperature T c . The second contribution ~ pbuik Lq originates from a microscopic 

drying layer near the solid boundary. At high temperatures the thickness Lq of the 
drying layer increases in accordance with the power law Lq ~ - In ( 1-T/T C ), indicating a 
drying transition at T c . The drying layer can be suppressed by strengthening the fluid- 
, surface interaction, by increasing the fluid-surface interaction range or by decreasing 

\^ ' the pore size. 
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1. Introduction 

Knowledge of the liquid density profiles near weakly attractive surfaces is necessarry 
for the understanding of various phenomena, such as hydrophobic attraction between 
extended surfaces in water pQ, slipping flow of liquids near a weakly attractive surface 

conformational stability of biomolecules in aqueous solutions, etc. Experimental 
studies evidence a depletion of the liquid density near weakly attractive substrates 
[31 IH El El 13 IE] • However, the shape of the liquid density profiles could not be derived 
unambigously from experimental data. 

In the vast majority of practically important situations, the liquid is close to 
equilibrium with the vapor, i.e. close to the liquid-vapor coexistence. At the liquid- 
vapor coexistence curve the vapor (liquid) phase undergoes a wetting (drying) transition 
near the solid boundary at some temperature T w (T^) E3 EEU H2j- Such a surface 
transition appears as the formation of a macroscopically thick liquid layer in the vapor 
phase near a strongly attractive surface (wetting transition) or as a vapor layer in the 
liquid phase near a weakly attractive surface (drying transition). When crossing T w 
(To) along the liquid- vapor coexistence curve, the liquid (vapor) layer could appear 
continuously (second-order or critical wetting/drying) or discontinuously (first order 
wetting/drying). 

Wetting transitions were extensively studied theoretically, experimentally and by 
computer simulations for various systems. The temperature of the wetting transition 
and its order were found to be strongly sensitive to the details of the fluid-fluid and 
fluid-surface interaction J2]- in particular, in model systems with short-range fluid- 
surface potential, the wetting transition can be of first or second order, depending 
on the strength of the fluid-surface potential [TQJ H31 Ej. The long-range fluid- 
surface interaction can change drastically the character and temperature of the wetting 
transition. Critical wetting transitions become first-order due to long-range fluid- 
surface interactions ^3]. The interplay between the short-range and long-range fluid- 
surface interaction potentials can suppress critical wetting up to the liquid-vapor critical 
temperature [TH] and produce a sequential (multiple) wetting transition [T7j . 

A drying transition which is accompanied by the formation of a liquid-vapor 
interface at some distance from the substrate was observed for short-range fluid-wall 
interactions (hard- wall [T%] IT§] |2*U] , square- well and truncated L J interactions 

[22]) in computer simulations and in density functional studies. Contradictory reports 
concerning the order of the drying transition in systems with short-range fluid-wall 
interaction were discussed in [23]. A long-range fluid-surface potential (e.g. via van der 
Waals forces) suppresses a drying transition at subcritical temperatures, which then can 
occur at the bulk critical temperature T c only [TSJ HHJ [22] ■ Since a long-range interaction 
between fluids and solids is unavoidable in real systems, the formation of a macroscopic 
vapor layer between the liquid and the surface is impossible. Indeed, a drying transition 
was never observed experimentally [2H] ■ 

Although a macroscopic vapor layer can not appear near a surface with weak long- 
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range attractive potential, an effect distantly related to a drying transition could be 
expected in the liquid phase below T c as the appearance of a drying layer near the 
surface (which is not a macroscopic vapor layer and could be considered as "embryo of 
drying" [2ED- This drying layer can noticeably influence the liquid density profile near 
a solid surface. 

Due to the necessary boundary conditions the behavior of the fluid density near a 
planar solid substrate can only be studied by computer simulations in a pore geometry. 
Most of the computer simulation studies of wetting and drying transitions were 
performed in NVT ensembles, where slitlike pores with symmetrical or asymmetrical 
walls were incompletely filled (see, for example Refs.[22[ I24j. In such simulations the 
average density of the confined fluid is deeply inside the two-phase region and a correct 
reproduction of the liquid-vapor coexistence is questionable. Even when the liquid- 
vapor interface in such a pore is well established, it is formed parallel to the wall and, 
therefore, the wetting phase is represented in such a system as a wetting layer only. 
In our computer simulation studies of surface transitions we use another approach: 
simulations of the liquid-vapor coexistence curve of fluids, confined in pores of various 
sizes, with subsequent extrapolation of the results to semi-infinite systems. 

The chemical potential of the liquid-vapor phase transition of the confined fluid 
is shifted with respect to the bulk. This suppresses the formation of a wetting or 
liquid layer along the liquid-vapor coexistence of the confined fluid below and above 
the temperature of the wetting transition, respectively [T2~] . To study the possible 
appearance of a drying (or vapor) layer in the saturated liquid near a weakly attractive 
surface, we have investigated the temperature evolution of the density profiles along the 
pore coexistence curve. The data obtained for fluids in pores of various sizes and with 
different fluid-surface potentials, are compared with the available theoretical predictions 
and experimental observations. 

A LJ fluid was confined in pores of width H = 40 a (a is the molecular diameter) 
with weakly attractive walls, which interact with the molecules of the fluid via a long- 
range potential. The well depth of the fluid-surface potential was about 70% of the 
well depth of the fluid-fluid potential. This is close to the interaction between neon 
atoms and a cesium substrate, the weakest known physisorption system, recently studied 
experimentally and by simulations [HJ EE] • In our recent paper |2Zj , the liquid- vapor 
coexistence curve of a LJ fluid was simulated in a smaller pore with H = 12 a) with 
the same fluid-surface porential. A depletion of the liquid density near the pore wall 
was observed in a wide temperature range along the pore coexistence curve without 
any trend towards the formation of a vapor or drying layer. It was found, that the 
fluid near the wall follows the universal power laws of the surface critical behaviour of 
Ising systems |2*%1 129j. In particular, a universal behavior of the local order parameter 
profiles, defined as the difference between the local densities p;(Az,r) and p v (Az,r) of 
the coexisting liquid and vapor phases at some distance Az from the pore walls, was 
observed in a wide temperature range |2Z1 EDI EH] • The intrusion of the surface critical 
behavior into the bulk fluid is governed by the bulk correlation length which diverges 
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Figure 1. Coexistence curve and diameter of the LJ fluid confined in the slitlike pore 
of width H — 40 a (circles). The bulk coexistence curve and diameter of the LJ fluid 
[27] are shown by solid lines. 

when approaching the critical point: 

£-(r) = eor-" (1) 

where r = (1 - T/T c ) is the reduced deviation of the temperature T from the bulk 
critical temperature T c , v = 0.63 is the universal critical exponent and £o is the 
system dependent amplitude. 

In view of the high sensitivity of the surface phase transitions to small changes 
of the chemical potential, the observation of a vapor or drying layer in a pore with H 
= 12 a was, probably, prevented by the shift of the liquid- vapor phase transition in 
such a rather small pore. To reduce the influence of confinement and to promote the 
observation of a vapor or drying layer, the pore width should be as large as possible. In 
the present paper we approach the bulk coexistence by performing simulations in a very 
large pore of width H = AO a. To extrapolate the results to a semi-infinite system the 
liquid-vapor coexistence of the LJ fluid in several pores of intermediate sizes between H 
= 12 and 40 a was simulated. 

2. Method 

Coexistence curves of a confined LJ fluid were determined using Monte Carlo simulations 
in the Gibbs ensemble (GEMC) GEMC simulations allow to achieve direct 

equilibration of the two coexisting phases, which are simulated at a given temperature 
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simultaneously in two simulation cells. The details of the simulations, parameters of 
the model fluid, as well as its bulk coexistence curve are given in our previous paper 
[27] . To study the behavior of the liquid near the surface, the L J fluid was confined in a 
slitlike pore with structureless walls. Each wall interacts with the particles of the fluid 
via the long-range potential of a single plane of LJ molecules: 

U w {z)=Aef [0A(a/z) 10 -(a/z) 4 }, (2) 

where z measures the distance to the wall and the parameter f determines the strength 
of the fluid-wall interaction relatively to fluid-fluid interaction and a is the diameter of 
the LJ particles. No truncation was applied to U w (z). In the present paper we report 
the results obtained for a pore of width H = 40 a and U w \(z) = U w (z, f = 0.3), which 
corresponds to a weakly attractive (solvophobic) surface with a well-depth of the fluid- 
wall interaction of about 70% of the well-depth of the fluid-fluid interaction. To study 
the effect of the pore size on the liquid density profile and the appearance of a drying 
layer, we also simulated the liquid-vapor coexistence and the liquid density profiles at 
T = 1.10 in pores of width H = 16, 22 and 26 a with the same fluid-wall interaction 
U w (z). 

To explore the effect of the strength of the fluid-wall potential on the liquid density 
profile we also simulated the liquid-vapor coexistence in a pore with slightly stronger 
fluid-surface interaction U w %(z) = U w (z, f = 0.4)- Additionally, the effect of the range 
of the fluid-wall interaction was explored by simulation of a system with a slower decay 
of the fluid- wall potential U w %(z), which follows from an integration of the individual 
interactions over the half-space of LJ molecules and can be described by the equation: 

U w3 (z) = 4e/* [(a*/zf - [a*/zf] , (3) 

The parameters of the potential U w s(z) were adjusted to get equal well depths of the 
potentials U w3 (z) and U w \(z). This was achieved using a* = 0.8328 a and f = 0.468 in 
equation (J3J). 

The density p used in the present paper is the reduced number density (scaled by 
a 3 ), while T is the reduced temperature (scaled by e/ks, where fcg is a Boltzmann's 
constant). The average density of the LJ fluid confined in a pore was calculated taking 
into account the volume accessible to the fluid molecules. The fluid-wall interaction is 
equal to zero at the distance 0.86 a from the pore wall and so in an operational approach 
this interval could be divided equally between the volumes of the fluid and the solid. 
As a result, the volume accessible to the fluid molecules is the total pore volume L 2 H 
reduced by a factor 1.022 for the pore of width H = 40 a. Accordingly, as the distance 
to the pore wall Az we used the distance to the parallel plane 0.5 a inside the pore. 

The density profiles of the liquid were obtained by Monte Carlo (MC) simulations 
in the NVT ensemble, using the average densities of the liquid phase obtained in the 
GEMC simulations. At high temperatures, the strong density gradient near the pore 
wall makes the determination of reliable density profiles in such large pores very time 
consuming. This problem was overcome by using two kinds of moves in MC simulations 
in the NVT ensemble. The first kind of moves is a standard MC move with a maximal 
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Figure 2. Temperature evolution of the liquid density profile along the pore 
coexistence curves of the LJ fluid confined in a slitlike pore with H — 40 a (upper 
panel) and H = 12 a (lower panel). 

displacement of molecule, which provides an acceptance probability of about 50 %. The 
second kind of move is a long-distance molecular transfer inside the simulation box: an 
attempt to place randomly chosen fluid molecule into randomly chosen position. This 
move is similar to the one used in GEMC simulations for molecular transfers between 
the two simulation boxes. Such long-distance molecular transfers essentially improve 
the sampling of density profiles, which show a strong gradient normal to the pore wall. 

3. Results 

The obtained densities of the coexisting vapor (p v ) and liquid (pi) phases and the 
diameter ((pi + p v )/2) of the coexistence curve of the LJ fluid in the slitlike pore with 
weakly attractive walls of width H= 40 a are shown in figure 1. The critical temperature 
of the pore coexistence curve is depressed due to the effect of the confinement. It 
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Figure 3. The density profiles of the coexisting liquid phase of the LJ fluids confined 
in the slitlike pores with H = 40 a (symbols) and 12 a (solid lines). The profiles for T 
= 1.04 and T = 1.12 are shifted to larger Az by 3 and 6 a, respectively. 

is estimated as T = 1.165 ± 0.005, that is only about 0.02 below the bulk critical 
temperature T c = 1.1876 [HHI2Zj- This indicates that in this case the pore coexistence 
curve is rather close to the coexistence curve of the bulk fluid. In such a wide pore, 
we may therefore expect that below the pore critical temperature the properties of the 
fluid near a wall are only slightly influenced by the presence of the opposite wall. 

Figure 1 shows that the average density of the liquid phase at the pore coexistence 
curve is noticeably lower than the density of the bulk LJ fluid at the same temperature. 
This effect becomes more pronounced at higher temperatures, i.e. at about T > 1.05 
and it originates from the depletion of the liquid density near a weakly attractive wall 
[2T| 130*1 EH HHH EH!- Indeed, a pronounced decrease of the liquid density near the wall 
is observed in the large (H = 40 a) as well as in the small (H = 12 a [22j) pore (see 
figure 2). At high temperatures the liquid density profile is flat in the pore interior 
(constant p;(Az,r) in the central part of the pore) and gradually decreases toward the 
pore wall, displaying only a single small oscillation near the pore wall (figure 2, upper 
panel). Upon cooling, the region of density depletion shrinks and additional density 
oscillations develop gradually on the density profiles. A qualitatively similar behavior 
of the density profile is seen in the liquid phase in the small pore of width H = 12 a 
(figure 2, lower panel). In the small pore, however, the density depletion spreads over 
the whole pore near the pore critical temperature and the flat part of liquid density 
profile is not observed. 

A direct comparison of the density profiles of the saturated liquid in the small 
and large pores is shown in figure 3 for some temperatures. At T — 0.80 the profiles 
are practically identical in both pores. The density oscillations are slightly larger in 
the small pore. At high temperatures, the density depletion is weaker in the small 
pore and the difference between p;(Az,r) in the small and the large pore becomes more 
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Figure 4. Scaling plot for the liquid density profiles in the small pore with H = 12 a 
for T= 1.04, 1.05, 1.06, 1.07, 1.08, 1.09, 1.10, 1.11, 1.12 and 1.13 (upper panel) and 
in the large pore with H = 40 a for T = 1.04, 1.05, 1.06, 1.07, 1.08, 1.09, 1.10, 1.11, 
1.12, 1.13, 1.14, 1.15 , 1.155 and 1.16 (lower panel). The highest studied temperature 
for each pore is shown by a thick solid line. 

pronounced with increasing temperature. 

It was shown [23 EI] that the local order parameter profiles Ap(Az,r) = (pi(Az,r) 
- p v (Az,r))/2, in the small pore of width H = 12 a with weakly attractive walls collapse 
into a single master curve, when the local order parameter is normalized by the bulk 
order parameter, and the distance to the pore wall is normalized by the bulk correlation 
length (see figure 9 in Ref. |27J). Such a universal behavior of Ap(Az,r) evidences 
that the drying transition, which is expected at the bulk critical temperature, does not 
influence the liquid density profiles in such a small pore. In the large pore of width H 
= 40 a, we do not obtain such universal behavior of Ap(Az,r) when using the same 
approach, as described above j2Z] • We have found, that this discrepancy originates from 
the quite different temperature dependences of the density profiles of the liquid phase 
in small and large pores. 

The liquid density profiles at various temperatures can be compared, using the 
normalized density pi(Az,r) / pi^uikij), where pi t buik(j) is the liquid density of the bulk 
LJ fluid (the values, obtained in Ref. [2Z1> are shown in the last column of Table 1). 
Because the intrusion of the surface perturbation into the bulk fluid is governed by the 
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Table 1. The values of the fitting parameters Lq and pi t buik in equation El when 
fitted to the liquid density profiles in the large pore of width H = 40 a. Uncertainties 
of the fitting parameters correspond to the confidence level 95%. The values of p® buU , 
were obtained by direct GEMC simulations of the liquid- vqpor equilibrium |27) . 
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bulk correlation length £, it is reasonable to display the distance to the pore wall in 
terms of £. Taking into account the universal temperature dependence of £ along the 
coexistence curve (equation (JTJ), we can use Az/r~ u as a normalized distance variable. 
The normalized density profiles of the liquid in the small and the large pore are compared 
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Figure 5. Liquid density profile in the large pore with H= 40 a at T= 1.10 (circles) 
and the fit to equation with the parameters shown in Tabic 1 (solid line) and to 
equation 0] (dotted line). Coordinates of the inflection points are shown by dashed 
lines. The inset gives a sectoral magnification. The value pi^uik found from the fit is 
shown by the horizontal line in the inset. 



in figure 4 in the temperature range from T = 1.04 up to the highest temperature of 
the two-phase coexistence in the pore. 

In the small pore, the normalized density profile pi(Az,r) can be roughly described 
by the exponential equation: 



The values of the parameter Iq obtained from the fits do not exceed a few tenth of 
a, whereas the value of the fitting parameter £ is about 0.3 a, close to the value of 
£o; obtained from the master curve of the order parameter profiles in Ref. j22j. In 
the large pore, the shape of the liquid density profiles is qualitatively different from 
that in the small pore (figure 4). Moreover, in the large pore the normalized profiles 
pi( Az, t) / pi^ u ik(r) vary with temperature non-monotonously with the "lowest" profile 
at T — 1.11 (see figure 4). The shape of the liquid density profiles at any temperature 
can not be fitted to the exponential equation (J3J). 

These peculiarities of the liquid density profiles in the large pore indicate the 
possible presence of some non-universal contribution to the surface behavior. As the 
increase of the pore size can be considered as an approach to a semi-infinite system, 
one may expect that this peculiar behavior is connected with the formation of a drying 
layer, caused by a drying transition at the bulk critical temperature ^2] • If such a drying 
layer appears in a liquid phase near the surface of a large pore, the liquid density profile 
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Figure 6. Temperature dependence of the thickness of the drying layer Lq and 
correlation length £_ , obtained from the fits of equation © to the part of the liquid 
density profile pi(Az,r), which extends from the pore center to various distances Az CMt 
from the surface: Az cut — a (squares), 1 a (circles), 2 a (triangles) and 3 a (stars). 
A logarithmic behavior described by equation Q with const = -0.05 a and -0.85 a is 
shown by upper and lower solid lines, respectively. The power law JD with £o = 0.28 
a is shown by a dashed line. 



should be described by an interface-like equation. The density profile at the intrinsic 
liquid-vapor interface, derived from the van der Waals theory, could be described by the 
following equation [37| 13%]: 

/ A \ Pl,bulk( r ) - Pv,butk{ T ) . u f Az — L \ Pl,bulk{r) + Pv,bulk( r ) / c s 

Pl {Az, t) = tank I — — I + , (5) 

where L is the distance from the liquid-vapor interface to the solid surface. This 
intrinsic interfacial profile can be affected by capillary waves, i.e. thermal fluctuations 
of the interface jHlj. Far from the critical point the bulk fluctuations are small and 
capillary waves are the dominant factor, which determines the thickness of the interface 
[4"U] . Upon rising the temperature the density fluctuations in the bulk phases become 
more important, yielding an intrinsic interfacial thickness proportional to the increasing 
bulk correlation length, as described by equation (jSJ). Since the expected interface is 
close to the wall, the capillary waves can be suppressed by a long-range fluid-surface 
interaction. 
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Figure 7. Variation of the fitting parameters pi,bulk and £_ at T — 1.15 (triangles) 
and T = 1.16 (circles), when equation 10} is fitted to those parts of the liquid density 
profile pi(Az,r), which extend from the pore center to various distances Az cut from the 
surface. The values of the liquid density at the bulk coexistence curve are shown by 
solid lines (upper panel) . The values of the bulk correlation length £_ , obtained from 
the fit of equation © to pi(Az,r) (Table 1) are shown by dashed lines (lower panel). 



We have found, that the liquid density profile in the large pore can be described by 
equation (JHJ at low temperatures, when the density of the bulk vapor is close to zero. 
However, equation (j3J) fails to fit the liquid density profiles at high temperatures, because 
the value of the density near the pore wall is even lower than the bulk vapor density 
at the same temperature. This result is consistent with the theoretical expectations, 
that a vapor layer never appears near a surface with long-range fluid-surface interaction 
[Tol Ho] . So, it is more reasonable to suggest not a vapor but a drying layer near the 
surface. Since the density distribution in a drying layer is unknown, we fitted the liquid 
density profiles to equation (JHJ) with p v ^uik as a free parameter p*. Such an approach 
provides a satisfactory description of the simulated profiles apart from the first layer of 
about er width, where a density oscillation, caused by fluid-surface interaction potential, 
is noticeable even at supercritical temperatures. The values of p* obtained from the fits 
were found close to zero for all studied temperatures. An example of such a fit with p* 
= for T = 1.10 is shown in figure 5. So, we propose the following equation to describe 
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the liquid density profile pi(Az,r) near a weakly attractive surface: 



(Az,r) = fi!f^[tan*(^!) + l 



(6) 



There are three parameters in equation (jfij): Lo, the distance of the inflection point of 
the such defined liquid-drying layer interface from the surface (thickness of the drying 
layer), the bulk correlation length at the liquid- vapor coexistence curve and pi^uikij), 
the density of the bulk liquid. The first density oscillation (which extends over about 
1 (j near the surface) is caused by the localization of the molecules in the well of the 
fluid-wall potential. Since it strongly deviates from the smooth density profile outside 
the first layer, it is reasonable to exclude the first layer from the fits. The values of the 
parameters obtained from the fits of equation (0) to the simulated liquid density profiles 
are shown in Table 1. 

The temperature dependences of the correlation length £_ and of the thickness of 
the drying layer L , obtained from these fits, are shown in figure 6 (solid circles). The 
correlation length shows a temperature dependence, which is rather close to the simple 
power law described by equation (JTJ) with an amplitude £o = 0.28 a (figure 6, lower 
panel, dashed line). Such an amplitude of the correlation length agrees well with the 
value £o = 0.3 a, obtained previously from the master curve of the order parameter of 
a LJ fluid in a narrow pore [2*7j . 

Another estimate of the amplitude of the correlation length can be obtained based 
on the exponential decay of the perturbation caused by the surface: far from the surface 
the deviation of the local density from the bulk value is governed by the bulk correlation 
length and should follow equation (jlj). Thus, the progressive exclusion of an increasing 
part of the density profile near the surface from the fit, should show a convergence of 
the fitting parameters, including The values of obtained from fits of equation 
(@J) to liquid density profiles, which extend from the pore center to various distances 
from the surface Az cut are shown in figure 7 (lower panel) for two temperatures. A clear 
convergence of the fitted values £_ is observed, when the fitted part of the density profile 
does not approach the pore wall closer than 6 a at T — 1.15 and 8 a at T — 1.16. These 
converged values of £_ are remarkably close to the values of obtained from the fit of 
the data to equation ©, shown by dashed lines (figure 7). 

The values pi^uik obtained from fits of pi(Az,r) to the interfacial equation (JBJ) (Table 
1) are almost equal to the liquid densities in the pore center, obtained by fits of the 
data near the pore center to the exponential equation (jlj) (see solid lines in figure 7). 
The values pi^uik obtained from these fits agree very well with the bulk liquid densities 
obtained by direct Gibbs ensimble MC simulations of the liquid- vapor equilibrium |27j . 
which are also presented in Table 1. 

Next, we test the sensitivity of the fitting parameters in the interfacial equation 
© from the portion of the liquid density profile used for the fits. For this purposes we 
fitted equation © to the part of pi(Az,r), which extends from the pore center to some 
distances Az cut from the surface. The fitting parameter pi^uik was found practically 
independent on the choice of the fitting interval. The values of the correlation length 
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obtained from the fits, depend on the choice of Az cut (from to 3 a) only at T < 
0.90 (see figure 6, lower panel) due to the strong density oscillations near the surface at 
low temperatures. At T > 0.90 the fitting results for the correlation length £_ are not 
sensitive to Az cut . 

The change of the thickness of the drying layer Lq with temperature is shown in 
the upper panel of figure 6. At the lowest temperatures (T < 0.90) the fitted value of 
L is determined by the choice of Az cut , namely L Q « Az cut . This correlation is clearly 
imposed by the strong density oscillations near the surface. Thus, if the thickness of a 
drying layer does not exceed 1 a, meaningful value of L can not be determined from 
the fits. At higher temperatures, L does not depend noticeably on the choice of Az cut . 
Two temperature regions could be distinguished in the temperature dependence of L Q 
at T > 0.90. In the temperature interval from T = 1.11 to 1.16 the thickness of a drying 
layer Lq clearly follows the logarithmic dependence: 

L = In (r _1 ) + const (7) 

with const « -0.05 o (see Fig. 6, upper panel, upper solid line). A quite similar 
logarithmic temperature dependence of L is observed in the low temperature interval 
from T = 0.85 to 1.06, but with another value const « -0.85 o. There is a crossover 
between two kinds of logarithmic behavior in the temperature interval from T = 1.07 
to 1.10. Note, that the observed logarithmic temperature dependence of L is similar 
to that predicted for the divergence of the thickness of the wetting layer in the case of 
a critical wetting transition [12J. 

So, equation |B1 allows the description of the liquid density profiles in a wide 
temperature range using two properties of the bulk fluid (pi t b u ik and £_) and the thickness 
of the drying layer L . The degree of the universality of the liquid density profiles can 
be illustrated by a master curve, which is essentially improved in comparison with figure 
4, where the presence of a drying layer was neglected. Now, we introduce a normalized 
length scale, which is measured from the liquid-drying layer interface , as (Az - Lo)/£_. 
Using the fitting parameters L and p^buik from Table 1, the liquid density profiles 
collapse on a single master curve in a wide temperature range 1.04 < T < T c (see figure 
8). The dashed line in figure 8 represents the interfacial equation (jHJ) in normalized 
coordinates. This means, that equation © allows a perfect description of the liquid 
density profiles in the whole temperature range, where the drying layer can be clearly 
detected. In the temperature range above T = 1.11, where the thickness of the drying 
layer follow equation (see figure 6), pi(Az,r) can be predicted if the bulk density, 
the bulk correlation length and the const in equation (jZJ is known. Note, that the latter 
value can be estimated from a single density profile. 

We also test the possibility to describe the density profiles of the liquid in the small 
pore with H = 12 a by the interfacial equation (JBJ). The obtained fitting values of 
Lq are about 1 o in the whole temperature range (figure 9). The fitting values of Lq 
in the small and large pores coincide at low temperatures (T < 0.90), where they are 
strongly influenced by the density oscillations and, therefore, are meaningless. As far 
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(*-L )/4_ 

Figure 8. Scaling plot of the liquid density profiles pi(A z,r) in the large pore H = 
40 a for the same temperatures as in figure 4. The parameters pijmik and Lq were 
obtained from fits of equation to density profiles pi(A. z,r), which were cut at the 
distance 1 a from the surface (Table 1). For T — 1.16 the value pi,bulk from the bulk 
coexistence curve was used. 

as a thickness of the drying layer below 1 o can not be detected from the fits, the real 
thickness of the drying layer can be essentially lower. Thus, the fitting values of Lq for 
the small pore with H — 12 a, shown in figure 9, can not be considered as an indication 
of a drying layer in this pore even near the pore critical temperature. Obviously, the 
fluid layer near the surface, highly localized in the well of a fluid-wall potential, can not 
be considered as a drying layer. This observation agrees with a previous analysis of the 
local order parameter in the same pore [2Z|, which indicated the absence of a drying 
layer. 

The values of £_ found from the fits of pi(Az,r) in the small pore with equation 
(JHJ) are noticeably smaller than £_ obtained from the fits of the order parameter in the 
same pore [21] and £_ in the large pore (figure 9). So, the liquid density profiles in 
the small pore can not be described correctly by equation (jSJ: the effective thickness of 
the drying layer is ultimately overestimated in such a pore, resulting consequently in an 
underestimation of the value of the correlation length. 

The drying layer, clearly seen in a large pore with H = 40 a, disappears in a small 
pore with H= 12 a and the same fluid-surface potential. This agrees with the theoretical 
expectation of strong sensitivity of the thickness of the drying layer to the chemical 
potential, which is shifted due to the confinement in the pore [T2J. One may expect, 
that the thickness of a drying layer Lq should continuously increase with increasing 
pore size and achieve some saturation value L Qjin f in a semi-infinite system. To check 
this expectation we have simulated liquid-vapor coexistence of the LJ fluid in several 
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Figure 9. Temperature dependence of the thickness of the drying layer Lq and of the 
correlation length obtained from fits of the liquid density profiles p;(Az,r) cut at 
the distance 1 a from the surface for large (solid circles) and small (open circles) pores. 
Solid and dot-dashed lines represents equation Q with £o = 0.28 a (this paper) and 
£o = 0.30 a [23 , respectively. Dashed line corresponds to equation (|7J with const = 
-0.05 a. 
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Figure 10. Liquid density profiles pi(Az,r) at T = 1.10 in pores of various width 
with identical fluid-surface interaction potential, described by equation @ . The liquid 
density of the bulk fluid is shown by a horizontal line. 
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Figure 11. Dependence of the thickness of the drying layer Lq and the correlation 
length obtained from fits of the liquid density profile p;(Az,r) with equation JBJ, 
on the inverse pore width 1/H (circles). Linear fits of the data are shown by lines. 

pores of various sizes with the same fluid-surface interaction potential at T — 1.10. The 
obtained liquid density profiles shown in figure 10 indicate, that the density depletion 
becomes more pronounced in larger pores. The fits of pi(Az,r) with equation © allow 
an estimation of the parameters L and £_ in all studied pores, and their dependences 
on the pore size are presented in figure 11. Note, that at this temperature, L exceeds la 
and the drying layer is indeed detectable in all pores except the smallest one with H= 12 
a. As the shift A/x of the chemical potential of the phase transition in a pore relatively 
to a semi-infinite system is proportional to 1/H |JT], it is reasonable to consider the 
dependence of Lq and on 1/H. Both parameters show essentially a linear dependence 
(see figure 11). Extrapolation of these linear dependences to semi-infinite system {1/H 
— > 0) gives the values: L 0jin f = (2.67 ± 0.15) o and £_j n / = (1-65 ± 0.04) a. Assuming 
that a drying layer effectively decreases the pore size to H - 2 L , the dependences 
shown in figure 11 become non-linear and their extrapolation to semi-infinite systems 
by second-order polynomial fit gives Lo,m/ = (2.79 ± 0.47) a and = (1-73 ± 0.09) 

a. The values of extrapolated to a semi-infinite system, correspond to an amplitude 
of the correlation length £ ~ 0.32-0.33 a. This value is about 15 % higher the value £o 
= 0.28 cr obtained in the pore with H = 40 a. 

Finally, we have tested the general theoretical expectation, that the formation 
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Figure 12. Upper panel: liquid density profiles pi(A z,r) at T = 1.10 (circles) in 
pores of width H — 40 a with U w (^(circles), U w \(z) (triangles) and V w i(z) (squares) 
fluid- wall interaction potentials. Lower panel: the same profiles shown as a distance 
to the inflection points (see text). 

of a drying layer could be suppressed by a strengthening fluid-surface interaction or 
by extending its attractive range. For this purpose we have simulated liquid-vapor 
coexistence and density profiles of the coexisting phases in the large pores (H = 40 
a) at T — 1.05, 1.10 and 1.15 with stronger interaction potential U m i(z) and with the 
slower decaying potential U w s(z) (see section Methods for more details). The obtained 
density profiles at T — 1.10 are compared in figure 12. Evidently, the drying layer is 
strongly influenced by the interaction potential. Namely, strengthening of the fluid- 
surface potential by about 33 % causes a shrinking of a drying layer at T — 1.10 from 
2.26 a to 1.76 cr, i.e. by about 20 %. At T = 1.05 the effective thicknes of the drying 
layer is suppressed to about 1 cr. An extension of the attractive range of the fluid- 
surface potential from ~ r~ 4 to ~ r~ 3 has an even stronger effect: the drying layer 
almost dissappears (Lq < 1.4 cr at T — 1.10). 

Starting from the inflection point, the shape of the liquid density profiles remains 
highly universal for all considered fluid-surface potentials, as it is determined mainly by 
the bulk correlation length (see lower panel in figure 12). Note, that the details of the 
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Figure 13. Temperature dependence of the excess desorption of the liquid in the large 
pore: total excess desorption T tot (solid circles) and its contributions due to drying 
layer (squares) and due to increase of the correlation length (triangles). Total 
desorption in the small pore is shown by dashed line. 



long-range attractive tail of a fluid-surface potential can produce some minor deviations 
of pi(Az,r) from the universal shape, described by equation JBj). An example of such 
deviations is also shown in the inset of figure 5. 

Using the obtained density profiles of the liquid, we analyzed the surface excess 
desorption, which describes a deficit of mass per unit surface area caused by depletion 
of the liquid density near the surface. The total excess desorption T tot contains two 
contributions: desorption due to the presence of a drying layer T L and desorption 
due to the density depletion Tg, governed by the correlation length. We calculated 
these two contributions by numerical integration of p/(Az,r) from the surface to the 
inflection point Tl and from the inflection point to the pore center T^. The total excess 
desorption T to t and its two contributions are shown in figure 13 as functions of the 
reduced temperature. For a semi-infinite system the excess desorption 1^ can be simply 
calculated by integration of equation El for the liquid-drying layer interface from the 
inflection point to infinity: 

r e (r) = pi >bu i k ln2. (8) 

This part of the total desorption behaves like ~ r~ v at any temperature and strongly 
diverges when approaching the critical point. Neglecting the density oscillations near 
the surface, the contribution can be obtained in a similar way using integration of 
equation El for the liquid-drying layer interface from the pore wall to the inflection point: 

r L (r) = PlMlk [L + 2£Jn (cosh(L /2^))} /2 (9) 

We analyzed the ability of equations (JHJ and © to describe the simulation results 
of the excess desorption. The excess desorption depends on the bulk fluid properties 
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Figure 14. Temperature dependence of the excess desorption of the liquid in double 
logarithmic scale: symbols are the same as in figure 14. Desorption and Tl described 
by equations [H] and and T to t — 1^ + are shown by dashed, dashed-dot and solid 
lines, respectively. 

only. Knowing the temperature behavior of the bulk liquid density pi^uik along the 
coexistence curve |27J and using the asymptotic equation (JTJ for the bulk correlation 
length with amplitude £o = 0.28 a, we can directly obtain the temperature dependence of 
from equation (jHJ) (dashed line in figure 14). Good agreement of the simulated values 
of with equation (jHJ) is observed in the whole studied temperature range. Deviations 
at low and high temperatures should be attributed to deviations of the correlation length 
from the asymptotic equation (JTJ) (see figure 6). 

The excess desorption T L , described by equation @, apart from the bulk 
parameters pi^uik and £_ contains also the parameter L , determined by the fluid-surface 
interaction. We obtain meaningful estimates of Lq for the temperatures, where the 
thickness of the drying layer L exceeds about 1 a. In this range L follows equation [7| 
(see figure 6) and therefore the contribution Tl varies roughly as ln(r~ 1 )(pi t b u ik)/2 , i.e. 
shows logarithmic divergence when approaching T c . The excess desorption Tl described 
by equation and the total excess desorption T tot = + Tl, described by equations 
E] and El are shown in figure 14 by dot-dashed and solid lines, respectively. Obviously, 
good agreement of the simulated values of T L with equation Q is observed in the range, 
where the temperature of Lq is known. 

In a wide temperature range the contribution from the drying layer to the excess 
desorption is dominant (figure 13). However, the contributions Tl and T^ approach each 
other with increasing temperature and practically coincide at r ^ 0.02 (figure 14). The 
contribution T^ obviously should dominate at higher temperatures, as the former 
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Figure 15. Temperature dependence of the thickness of the drying layer Lq (solid 
circles) and correlation length £_ (open circles) , obtained from fits of the liquid density 
profile p; (A z,r), which extend from the pore center to the distance 1 a from the surface. 
The logarithmic law described by equation Q with const = -0.05 a is shown by dashed 
line. The power law Q with £o = 0.28 a is shown by solid line. The temperature, 
where Lq w £_ is marked by a star. 



diverges as r whereas the latter diverges logarithmically only. 
4. Discussion 

We have studied the temperature evolution of the density profiles of a liquid along the 
liquid-vapor coexistence curve near a weakly attractive surface. To approach the fluid 
behavior in a semi-infinite system, the liquid-vapor coexistence curve of a LJ fluid was 
simulated in an extremely large slitlike pore with H = 40 a. Additionally, at some 
temperatures we simulated the liquid- vapor coexistence in pores of various sizes. Strong 
depletion of the liquid density near a weakly-attractive pore wall is observed. The shape 
of the liquid density profiles evidences the formation of a drying (not a vapor) layer near 
the pore wall with increasing temperature. This drying layer can be detected at T > 
1.0 ( r < 0.15), when the distance between the inflection point of the density profile 
and the pore wall exceeds at least one molecular diameter. We propose to describe the 
liquid density profiles near weakly attractive surfaces by equation (0), which assumes 
an interface between the liquid phase and the drying layer and allows estimation of the 
bulk correlation length £_ and the thickness of the drying layer L . 

The fits of pi(Az,r) with the interfacial-like equation (JHJ) and with the exponential 
equation (j3|) which describes the decay to the bulk density far from the surface, give 
similar values of the amplitude of the correlation length £o ~ 0.28 a. This evidences a 
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self-consistency of the proposed treatment of the liquid density profile and also indicates 
a negligible effect of the capillary waves on the studied interface between the liquid and 
the drying layer. This conclusion corroborates general theoretical arguments, concerning 
the influence of a long-range fluid-surface potential on the interface between vapor and 
wetting layer [T2j. 

The thickness of the drying layer L follows a logarithmic temperature dependence 
(J2J), starting from T w 1.1 (i.e., at r < 0.07). Such a logarithmic divergence with 
approaching wetting (drying) temperature is expected for critical wetting (drying) [T2j. 
This result is consistent with the existence of a second-order drying transition at the bulk 
critical temperature for long-range fluid-surface interactions jTH] and with experimental 
evidences for the absence of a drying transition at subcritical temperatures even in a 
system with extremely weak fluid-surface interaction jSlEH], similar to the one, used in 
our studies. Additionally, we have shown how the thickness Lq of the drying layer can 
be strongly suppressed by strengthening of the fluid-surface interaction, extension of its 
attractive range and by confinement. 

In the temperature range of our simulation studies, the correlation length does 
not exceed the thickness of the drying layer Lq (see Table 1 and figure 15). Since 
the divergence of the correlation length £_ is much stronger than the logarithmic 
divergence of L , when approaching the critical point £_ should exceed L Q . At the 
characteristic temperature t fh 0.01, shown by a star in figure 8, ~ Lq. Above 
this temperature the interface between liquid and drying layer approaches the surface 
in terms of correlation length. The effective "thickness" of the total region of density 
depletion near the surface should diverge asymptotically as the correlation length. This 
corroborates the theoretical expectations for the thickness of the wetting (drying) layer, 
when the critical wetting (drying) transition occurs at the bulk critical temperature 
|16j . Experiments with fluid mixtures indicate, that the partial wetting layer decays 
exponentially with distance from the surface and its effective "thickness" is proportional 
to the bulk correlation length £_ S3] • 

We found that the total excess desorption due to the depletion of the liquid density 
near a weakly attractive surface should diverge as the bulk correlation length ~ t~ u 
— t~ ' 63 when approaching the critical temperature. This strong divergence, described 
by equation (JHJ), originates from that part of the liquid density profile which extends 
from the bulk liquid to the inflection point of the liquid-drying layer interface. Slight 
variations of the fluid-surface potential do not affect noticeably this part of the profiles 
(see figure 12, lower panel) and, accordingly, the value of T^, whereas the thickness of 
the drying layer is highly sensitive to the details of the fluid-surface potential (figure 12, 
upper panel). 

The observed divergence of the excess desorption is stronger than V ~ t^~ v = 
r -o.305^ eX p ec ted from the theory of critical adsorption [13]. We do not see the possibility 
of a crossover of F to the power law of critical adsorption with further approaching 
the critical temperature, as the inflection point of the liquid-drying layer interface 
asymptotically approaches the value p c /2. Such a crossover could be possible only if the 
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drying layer transforms into a vapor layer and, accordingly, the density at the inflection 
point of the liquid- vapor interface approaches p c with increasing temperature. However, 
this scenario is forbidden for systems with long-range fluid-surface interaction, where 
vapor layer do not appears below the critical temperature [13 HH1 121] • Note, that recent 
experimental studies of the adsorption in binary mixtures indicate, that at subcritical 
temperatures T diverges stronger than ~ r^~ v |43| H3] . Such behavior may have the 
same origin, as that observed in our computer simulations. 
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